A&A manuscript no. 

(will be inserted by hand later ] 



Your thesaurus codes are: 

missing; you have not inserted them 



ASTRONOMY 

AND 
ASTROPHYSICS 



Theoretical interpretation of the X-ray properties of 
GRB960720 



F. Daigne and R. Mochkovitch 

Institut d'Astrophysique de Paris, CNRS, 



On 

9^ ' Received **.**.** 

0\ ■ 



bis boulevard Arago, 75014 Paris France 



accepted 



in 



> 

CO 

o 

On 
On 

Of 

6 

C3 



Abstract. BeppoSAX observations of the single pulse 
burst GRB960720 have allowed a detailed study of its X- 
ray properties: pulse width in different energy bands, spec- 
tral evolution from 2 to 700 keV, etc. We show that the 
early (0— 5s) X-ray emission and the gamma-rays are well 
explained by internal shocks in a relativistic wind while 
the late (5-20s) X-ray emission could come from the re- 
verse shock generated in the wind when it interacts with 
the external medium. The results for a medium of uniform 
density are compared to the observations. 
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1. Introduction 

The gamma-ray emission from GRBs is probably pro- 
duced by internal shocks in a relativistic wind whereas 
the afterglow (from X-rays to radio bands) is due to the 
external shock, i.e. the forward shock propagating in the 
ISM beca use of its interaction with t he wind (Rees and 
Meszaros, 1994 ; Wijers et al., 1997). Simultaneously, a 



reverse shock propagates in the wind itself. We illustrate 
here the possible contribution of this reverse shock to an 
X-ray emission perduring immediately after the gamma- 
rays. Such an emission has for been observed in the first 
GRB detected by Beppo-SAX, GRB960720. We use the 
detailed X-ray data made available by SAX for this burst 



(Piro et al., 1998) to make a comparison with our theo- 
retical results. 



2. The internal shocks 

GRB960720 has been observed both by BATSE and 
Beppo-SAX. It is a single-pulse burst, with a "FRED" 
profile. Its duration in the 50-700 keV band is around 2-3 
s but the X-ray emission lasts longer: Piro et al 



( P-998D 

show that the power-law between the width of the pulse 
and the energy (already known in the gamma-ray range) 
is observed down to 2 keV. They find W(E) oc E~ aA6 . 
We use a simple model to simulate internal shocks and 
build synthetic bursts: all pressure waves are neglected so 



that we consider only direct collisions between solid lay- 
ers. In the shocked material, the magnetic field reaches 
equipartition values (10 - 1000 G) and the Lorentz factor 
of the electrons is obtained from the dissipated energy per 
proton e using the formula given by Bykov and Meszaros 
( |1996| ) who suppose that only a fraction £ of the electrons 
is accelerated: 



C, m e c 2 



(a M ~ 0.1 -> 1 and 1 < a < 1.5) (1) 



For £ ~ 1 the usual equipartition assumption yields values 
of T e of a few hundreds: the gamma-ray emission is due 
to inverse Compton scattering on the synchrotron pho- 
tons. Smaller values for the fraction of accelerated elec- 
trons (C < 10 -2 ) lead to larger Lorenz factors (r e of a 
few thousands) so that the gamma-rays are directly pro- 
duced by the synchrotron process, which is the assumption 
made here. Internal shocks have been shown to success- 
fully reproduce the main temporal and spectral properties 
of GRBs (Daigne and Mochkovitch, |199S| ). 
We model GRB960720 with a wind emitted during 4 s 
and consisting in a slow and a rapid part of equal mass 
(see figure [j]). Two internal shocks are generated and we 
sum both contributions to the emission to construct the 
synthetic burst. The profile in the SAX 50-700 keV band 
looks very similar to GRB960720 as can be seen in figure 
0. However the X-ray emission does not last long enough 




Fig. 1. 

wind. 



Initial distribution of the Lorentz factor in the 
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so that the power-law relating W(E) and E is not repro- 
duced in this spectral range. 

3. Effect of a medium of uniform density 

We now consider the effect of the ISM, whose density n 
is supposed to be uniform. The external shock (forward 
shock) produces the afterglow. Simultaneously a reverse 
shock propagates into the wind. Its strength is comparable 
to that of the internal shocks, which are midly relativis- 
tic while the external shock is initially very strong and 
relativistic. We therefore adopt the same assumptions to 
compute the emission of the shocked material behind the 
reverse and the internal shocks. 

It is possible to derive a critical ratio of the total energy 
injected over the density (see Sari and Piran, 1995) : 



— m„c 5 f 8 T 3 



(2) 



for which the reverse shock will interfere with the internal 
shocks. Injecting the typical values used in our example 
(Lorentz factor f ~ 300 and duration T = 4s) gives 



E mj 
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(3) 



where E^ 3 



tfsttt — ; — • Assuming an efficiency f~ for 

10^/4^ erg/sr & > 

the conversion of wind energy into g^amma-rays by inter- 
nal shocks allows to obtain E %n "* — -f- from the observed 
burst energy E^. 

The total energy being fixed, if the density is smaller than 
the critical value, the reverse shock does not contribute in 
gamma-rays and only produces a delayed X-ray emission 
with an intensity which increases with n. If the density 
reaches or exceeds the critical value, the X-ray and the 
gamma-ray profiles become very affected. 
We have represented in figure ^ the profiles for a ratio 
E52 l n ~ 0-5 ( n ~ 0.15n c ,it). The gamma-ray profile 
is unchanged and the X-ray profiles are improved. The 
corresponding E - W(E) diagram in figure || shows a 
well reproduced power-law over the complete energy range 
with an index —0.45 in agreement with the observations 
of GRB 960720. In the spectrum the contribution of the 
reverse shock to the late emission appears like a X-ray 
plateau, which is observed (and is even more extended) in 
GRB960720. 



is reproduced) and the spectrum. 

It is now important to study the case of more com- 
plex environments. In hypernova models, the progenitor 
is a massive Wolf-Rayet star with a dense wind (M ~ 
3 10~ 5 M & /yr and Voo ~ 2000 km/s typically) leading 

to high densities (n ~ -. — ) near the source. Our 

first calculations show that in this case even the gamma- 
ray profile could be strongly affected by the reverse shock, 
which may represent a potential problem. 

References 

Bykov A. and Meszaros P., 1996, ApJ 461, L37 
Daigne F. and Mochkovitch R., 1998, MNRAS 296, 275 
Piro L. et al., 1998, A&A 329, 906 
Rees M.J. and Meszaros P., 1994, ApJ 430, L93 
Sari R. and Piran T., 1995, ApJ 455, L143 
Wijers R.A.M.J., Rees M.J. and Meszaros P., 1997, MNRAS 
288, L51 



4. Conclusions 

The reverse shock propagating into the wind following its 
interaction with the ISM has been shown to produce, for 
a sufficiently dense medium, a late X-ray emission. In the 
case of GRB960720, whose gamma-ray properties are well 
explained by internal shocks, taking into account this ad- 
ditional contribution greatly improves the X-ray profiles 
(the power-law between the profile width and the energy 
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Fig. 2. Results obtained for a ratio E^ J jn = 0.5. For E 1 = 10 51 /47r erg/sr and / 7 = 0.05 it corresponds to 
n = 4 cm~ 3 . Left panel : Profiles of the synthetic burst in the X- and gamma-ray bands. The full line takes into 
account the presence of the ISM whereas the dotted line only shows the contribution of the internal shocks. Middle 
panel : Pulse width as a function of energy. The dotted line corresponds to W(E) oc E~ 0A5 . Right pannel : Burst 
spectrum. The full line takes into account the ISM and the dashed line does not. 



